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{- ‘ ABSTRACT 


The effects of a neutron charge distribution are considered. 


It is suggested that the charge distribution is of the form 
Ex - 
( ee oo) i. 
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The effects of this charge distribution for various values 
of the parameters c and f, on the mean square radius of carbon, 
and on the charge form factor of carbon at high q are shown. The 
effect of the neutron becomes appreciable in carbon above the 
diffraction minimum. The neutron in this range of q” adds approx- 


imately 50% to the carbon form factor. 
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1. INTRODUCTION 


A review of current literature concerning the nuclear charge 
distribution shows that the effect of the neutron is largely ignored 
except in deuteron studies. This paper proposes to suggest what 
effect, if any, a finite neutron will have on the distribution of 
charge in any nucleus. 

Before the charge distribution in a nucleus can be calculated, 
the wave functions describing the motion of the nucleons and the 
nucleons' internal charge distributions must be known. The proton 
charge and magnetic moment distributions and the neutron magnetic 
moment distribution are known. The neutron charge distribution is 
not well known, however. 

If wave functions are known, then the neutron charge distribution 
may be related to the measured nuclear charge distribution. In this 
paper the diffraction minimum observed in electron scattering by a 


cle is used to determine the nucleon wave functions. 


2. THE NUCLEAR FORM FACTOR 


In order to develop a theory for the charge distribution 
of a nucleus, it is assumed that each nucleon maintains its 
identity in the nucleus. This assumption leads to models which 
successfully predict various experimental results. The assumption 
implies that each nucleon in the nucleus may be described independ- 
ently of the other nucleons. The nuclear wave functions may be 
written as follows: 

PR) POR) oo £1) 

(7) 
ale 
Al 
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i. describes the position of the center of mass of the 
ee particle. The W' implicitly include a spatial wave function, 
hereafter referred to as), Spin and isospin functions. All of 
these functions are orthonormal. The use of the othonommality 
Will be implicit in what follows. The total wave function, Y F 
is also normalized and is antisymmetric on interchange of any two 
nucleons (a consequence of the nucleons being fermions). 

The purpose of this paper is to examine the nuclear charge 
distribution @(R). This cannot be measured directly. What can 


be measured is the differential cross section for scattering electrons, 


ag 
re 


The cross section may be separated into a term due to the 
charge distribution, and a term due to the magnetic moment of the 
nucleus [1]. The charge distribution term leads to the charge form 
fagtor, F (q?) , Which, in Born approximation, is the Fourier trans- 


form of the normalized charge distribution. 


> = 
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The Born approximation will be used throughout this paper. 


The nuclear charge distribution = (R) , 1s determined by the 


nuclear wave functions and the charge distribution of the nucleus. 
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where if is the nucleon charge operator: 


(2~4) ? bs < goek-) | 


ee 
The f. are the charge operators for the individual nucleons. 


Putting equations (2-3) and (2-4) into equation (2-2) yields 
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Let Y. ce) be denoted by a Note that after integrating over all 


a coordinate, we have 
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We see from equation (2-5) and (2-6) that, if the single 
particle wave functions are orthogonal, the charge distributions 
for all of the protons are equal, those of all the neutrons are 
equal and the charge distributions are additive. 


Using (2-6) in (2-5) implies 


(2-7) F(9") = 


I 
ry |- 
== 
ran 
i 
we 
~ 
_ 
ma 
l- 
bo 
eo N 
oe 
we 
— 
—— 
ae 
3 
Pos 
ce — 
ou 
<-) 
Ze 


The relation between R andr. isR=7. +Yr.'. 


Insert the following equation into (2-7): 
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Equation (2-7) becomes 
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Integrate over R, 
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Split the sum into one over protons, f. = a asin one over neutrons, 


f. = f 
a n 
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Since we have integrated over ¥, Fis = F, for all i. Use (2-12) 


and (2-13) in equation (2-11). 
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Then write equation (2-14) as 
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The nuclear form factor has been separated into three parts: 


(1) Ey! the form factor of the proton 
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(2) Fue the form factor of the neutron, and 
(3) E areo a form factor due to the distribution of the centers of 
mass of the nucleons. 
F is valid for any nucleon distribution, id (R); i fe (R) 
cm-A 
may be expanded (since any function of angle may be expanded in terms 


of spherical harmonics) as follows: 
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where ye is lin (R) averaged over all angles or the spherically 
symmetric part of /P®) ; 

F(q*) may be expanded in terms of the multipole expansion of the 
charge distribution. This yields an electric monopole term plus higher 
order multipole terms in the form factor. The simplest nuclei to 
analyze are then the spherically symmetric nuclei, since these have 
only the electric monopole term. All further work in this paper will 


deal only with this type of nucleus. 
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3. MEAN SQUARE RADIUS DETERMINATION 


A quantity which is experimentally easy to determine is the 
mean square radius. We show below that, theoretically, this radius 
may be expressed in terms of the particle wave functions and in terms 
of the charge distributions of these particles. 

The form factor for a spherically symmetric charge distribution 


3(r) is 
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We have on expanding sin (qr): 
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From equation (2-16) we have 
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Using this in equation (3-3) yields 
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In equation (3-4) we have used 
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The directly measurable quantities in equation (3-4) are 


» 


the mean square radii of the nucleus and of the proton. The mean 
square center of mass radius is dependent on A and the configuration 
of the nucleons in the nucleus. This configuration may be different 
for different elements of equal A. 

It is thus possible to write the nuclear form factor in terms 
of the nucleon form factors and nucleon wave functions which describe 
the states of the nucleons. See equation (2-16). However, starting 
with experimental data on the nuclear form factor, it is impossible 
to separate the contributions due to the nucleon charge distributions 
and that due to the nucleon motion unless a model for the nuclear or 
nucleon structure is assumed. That is, the analysis will be model 
dependent. The validity of any model is measured by how well it 


agrees with experiment. 
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4. POSSIBLE NUCLEON CHARGE DISTRIBUTIONS 


There has been a great deal of work done on the proton 


which suggests that it has a gaussian charge distribution [2]. 
2 


> 
wy £0) 2 es 


at” rid 
The form factor corresponding to this charge distribution is 


id 
i + 
(4~2) fej. - € 


Figure 1 shows a comparison of this form factor to experimental 
data [2]. We assume here that f 1s correct, and that the charge 
distribution and the matter distribution are identical. With this 
assumption it seems valid that the matter distribution in the neutron 
could also be distributed in the same manner. That is to say, the 
matter distribution of the neutron is a gaussian with approximately 
the same radius as the proton. It is known that the total charge 
of the neutron is zero. However, itS magnetic moment suggests that 
the neutron does have a charge distribution. By analogy with the 


proton the following form is appealing: 
a 
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This is a double gaussian, each with unit charge. The negative charge 


has the larger radius. 
Our assumption that the matter distribution of the proton is the 


same as that of the neutron implies that the radius of the neutron's 


negative (larger) gausSian is approximately that of the proton. 
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We take f = .44F (slightly larger than the proton). This leaves 
one parameter to fix in the neutron. The free parameter will be 
fixed in this paper by using the harmonoc oscillator model and 
experimental carbon data. (For alternative method see Appendix C.) 


The form factor correspond to f is, in Born approximation 


5 ; a t 


Figure 2 is a for various values of o compared to experimental 
data [3]. This data does not deny any of the suggested curves, but 


favors F = 0. 
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Figure 2 
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5. HARMONIC OSCILLATOR WAVE FUNCTIONS 


In order to determine< re we are forced to assure a model 
to predict <x“ 7 . The harmonic oscillator potential is a simple 
potential which appears to agree well with experiment. This theory 
has the added advantage that its parameters are fixed independently 
of the mean square radius of the charge distribution. 

Consider non-interacting particles in a spherically symmetric 
harmonic oscillator potential well: 

| Us 

(-1) Ver) > _ ve 
The stationary states or orbitals occupied by the particles are 


solutions of Schrodinger's equation. 
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Using equation (5-1), and defining a =e equation (5-2) can be 


Wrileeen as 
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Using equations (5-4) and (5-5) in equation (5-3) 
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The solutions to the angular equation are the spherical harmonics. 


(See Appendix B.) 


For convenience in solving the radial equation, define 


(5-7) Riv) = a 


The radial portion of equation (5-6) may be written as 
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The solution of equation (5-8) and equation (5-7) yields: (See 


eax B.) 2 : 
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These equations also yield: 
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Where /\ = 2N + 1 denotes the energy levels. The harmonic wave 


functions are then (See Appendix B.) 
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The above solutionsare based on the central potential having 
a fixed origin. In reality, the central potential should be measured 
from the mean position of the nucleons. it has been shown that this 
complication introduces a factor € Ta into the form factor [4]. 


From equation (2-4) and the above consideration we have (See 


Appendix B.) 
’ é 
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From equation (2-16) we have 
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If the harmonic oscillator model is applicable and if the 


nucleon form factors F and Ey, and the nuclear configuration (1.e., 


n, 1, m must be known for each nucleon) are known, then the total 


nuclear form factor may be written. For example, for p-shell 


nuclei (4¢A$16, 2¢Z<8) [5]: 
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from Appendix B have been used. In Carbon A=12, Z=N= 


tt 


A 
The Die a 


a Ot ay 
igi \ = th (lm 
eS. ore 


ie 


Or the total form factor is tt is sa ( -+ ) 
— 3, 
Bcty = (FR Cie tae 


2 has a zero at q” = * . Therefore, our cross section measurements 
a 


should also show a zero at this point. Experimental data show a 


os 


Marked minimum in the cross section which can be used to fix a ‘ 


th 
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6. COMPARISON TO EXPERIMENTAL DATA 


In the harmonic oscillator model, the center of mass form 


factor, Eerie is completely determined by the position of the 


minimum in a plot of Be (q°) vs. qi”. 
Carbon offers a very good method of fixing the neutron parameters 
since it is spherically symmetric and the experimental data show a 


a 


prominent diffraction minimum [6]. 
=—_ -) 
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Figure 3 is a comparison of equation (6-1) to experimental values 
for the neutron parameters of Figure 2 [6]. This figure indicates 
that the best neutron form factor is given by 
; $63) a C44) | 

+ ~ e. é| 


(6-3) » = ¢ 


The mean square radius of carbon is then 
2 1 2 49." % 
a 2 < ‘. P+ a”, 24 : 


(6-4) 1 a 2 Q 
1A ee! F eee GHEE 


This radius is, of course, model-dependent because of the manner 
in which&r® 2. has been found. However, it agrees well with current 
model independent determinations of < 17 c at Darmstadt in which 


fr°Ac = .587 + .19 F* [7]. 
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7.) DISCUSSION 


The preceeding results suggest two altermative approaches 
for determining the validity of the neutron radius in equation 
(6-4). First, assume that the nucleons are-described by harmonic 
oscillator wave functions, in which the range parameter a is adjusted 
to give the experimentally observed diffraction minimum in F (q”) : 


This allows us to write 
. 

— N a6, ee ee 

We would then use the model independent <r*/ determined at low gq? ; 

the accepted value of < ms 7 determined from Ht data, and the shell 

model prediction of < re. . Extreme accuracy in the measured 

quanti is required due to the small size of< > : 

The second alternative would be to assume that the nucleon 
center of mass distribution is the same for elements with equal A 
but different Z. This assumption is based on the charge independence 
of nuclear forces. The assumption could be checked in a model- 
dependent way by seeing whether or not their diffraction minima 
occurred at the same value of q. Measuring the model- independent 
radii of two elements with the same atomic weight at low g and 


subtracting the following: : . 
evR = < ‘? e LP + Sema? 


= | ra *¥ r 
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Thus< r“7 is obtained in terms of only the model independent 
‘mean square radii. The assumption, however, that< aa =< mee 
would only be valid if the nucleon configurations of the elements 
were identical. One possibility would be to use chromium 50 and 
titanium-50. These two nuclei have total angular momentum (J) equal 
to zero. Therefore, they have only a monopole form factor. Any 
other pair of mirror nuclei would do equally well. 

Both of these methods could be checked by comparing the 
theoretical form factors to the experimental form factors at high 
qi’. This is the manner that the neutron parameters were fixed in 
this paper. | 

The limitation on the parameter f in the neutron charge 
distribution (that the neutron matter distribution equals the 
negative charge distribution) is of questionable validity. If 
this restriction is removed, and we use instead the electron- 
neutron interaction condition, 

(7-5) oF, =e -} oe ol [8], [9], 
) 4! Q=0 

which implies 
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Then the difference between the two neutron parameters is .084 FO) 
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When this restriction is used and the form factor 
gta 7. 
S - C by = . iy 1s again fitted to various 
experimental data, a good fit is obtained for 6° = .09 Fe 
f° = 174 F* . 

This choice of parameters leads to a form factor which gives 
a smaller radius at low a’, but at a of about 25 pre the radius 
is much greater than with the previous choice of parameters. The 
form factor of the neutron, with this choice of parameters, becomes 
about three times that of the proton at qi” = 25 a This doesn't 


seem reasonable. The previous selection of neutron parameters is, 


consequently, preferable. 


a 


8. CONCLUSION 


Current experimental data point to the possibility that the 
charge distribution of the neutron has a greater mean square radius 
than previously thought. It appears that a value es ge =, 25° 
is not unlikely. In fact, such a value is supported by current 
measurements of nuclear radii. 

It should be noted that it is presently impossible to separate 
the effect of the neutron from the effect of the wave functions in 
currently available experimental data on electron scattering. A 
modification of the wave functions can give form factors which fit 
the data over a large range of q without considering any neutron 
effect. 

The approach suggested in this paper has the advantage of 
simplicity. Simple wave functions and a simple neutron give agree- 
ment with PG’) data. 

If the charge distribution of the neutron can be determined 
independently (for example, through neutron-electron interactions) , 
then this information could be used in fitting wave functions to 


experimental data. 
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HARMONIC OSCILLATOR WAVE FUNCTIONS AND FORM FACTORS 


Spherical Harmonics 
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Electric Monopole Form Factors for Independent Particles 
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APPENDIX C 
CLASSICAL MAGNETIC MOMENT DETERMINATION 
Classically the magnetic moment m is defined as 


—>» _*. 
= on x 
(c-1) M 93 we RU Ce 


one [10} 

where Y is the distance from some point of interest and J is the 

current density. For an arbitrary charge distribution, id (x) , rotating 
NN 

about an axis with velocity dJ) « 
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Using this in equation (C-1) yields 
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this becomes 
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We are interested in finding the magnetic moment for neutrons 


and protons projected onto the R axis. For a spherically symmetric 


charge distribution we have 
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The second term integrates to zero when the angular integration is 


per formed. 
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So the magnetic moment about the Z axis is for the proton 


Co an Ae 
r. 

and for the neutron 

(C30) Mie wie mS gr? 
2 é. ‘ 


The ratio of the proton moment to the neutron moment is then 


- r 
_ we @, lo a a 
Ne wv b, \ a 


assuming tJ, = Jn 


This agrees in sign with the value predicted ford 7 as 
compared to< x 7 but is only the correct order of magnitude and, 
therefore, gives neutron parameters larger than those in this 


paper. 
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